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We indicate how to extend the 3iisch slgorithn to handle £ class of 
special functions defined in ten*s of integrals. Host of the 
integration aachinery for this class of function* is siaiiar to the 
:^achir.ery in the algorithm which handles io^arithns. A program 
embodying much of the extended integration slgoriti^-m has been 
written • Ii u&s used to check & table of integrals and it succeeded 
in finding sona misprints in it. 



Tnc Integration of a Class of S?ccicl puncciMd flu Ghs KU*ft .U t o:i:;„.. 



The purpose of this mcaorcnCur* ia to >cinc out the: similarity for the 
purpose* of integration of tha logarithm fu&ctioa and a cU&s o£ special 
functions defined ia teraa of ta»sr*La. Xccc LtfAG 



log(h(x)) -J -de. 



that is, tha logarithm can be defined &s an integral. Kote also that we do 
not restrict ourselves to a single function log<x)j but consider the family 
of functions logC^(^))- Here we shall consider the class of spacial func- 
tions F defined by 
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<h(s)) -J g<c)dt, 



where £ is a constant, and g<t) can be obtained by rational operations upon the 
variable z and by exponential and logarithmic terzis involving t. The function 
£ ;r*3y also contain special function terms not involving F, but it cust not 
5a i^tegrable in terns of the allowed functions. This class of special 
functions inducts 
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the error function errfhfc;)) ° t" ! e ct, 

the exponential integral 2j\h(x)) ■ 7— dt, 
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ane. ^.c :>pancc . ta «t — *x. 
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Ve clain chat ituch of e.io ;r-ichiniry icr i^tc^r^tin^; cruras* Il-v; &nva*¥ittii 
such specie! functions is GMBinod in the Io;;aric^c. cmc of the Itisch 
algorithm fi]» An extension to civc hisch ai^orich^ is to bti preferred to 
traditional integration motbods (eg* , integration by parts, substitution 
of variables) because it is *. decision procedure. ?hts a-enns iu: the 
procedure should be Able co orcein tie j&ugnl if in is expressible in 
closed form or to show thiC the integral ac&MG be so ezvpressed. VradieioMi 
xerhecs cannot oekt such £ decision and are particularly weslt if the 
InzegKcad involves division in a noatrivicl wy« or. the other hand thfc 
Xisch al^orithsi we shall employ i$ WO ak in ISiEt ic does cot handle algebraic 
expressions (e.g., /*0* However, such a restriction hes been rexoved by 
S£sch in extensions to cha algorithm described in [3»4J. 

Our description of the seeps necessary co extend the Risch ai&oritha 
«ill concentrate on the particular case of the error function. We have 
extended the current sute of our implementation of the RUch algorithm co 
habile the error function. At the present time cur tapl*man«tion of the 
2ci£oh cigCEithni contains a nearly complete implementation of the logarithm 
ca-re and a partial icsplenantation of the exponential caae* Details about the 
appsawfc t£.<ea in our implementation ire to bs cc^nd in [2 J. 

the fine thing one notices in att3=?ting an extension is that the 
introduction of special functions will cause a modification to be Bftde M 
Ciw Licuvillc theorem which is basic to the 2£sch algorithm, the Liouville 
thecren says i roughly speal;ing» that the ir.ce^ral of an expression 5 i^ the 
-w c£ in expression whioh is rational in the logarithmic and exponential 
tirr*s already in 2 and constant t:uLcipLas of new io^&rizh-^ic certis. "he 



thar ciiJarbariarica ^ooj cat wksmm , 4t _ ta ;; < »*iyifci&l or Ic^arCi^ic 
tersu,, y-oriov^r, diiicrcaz;_atiG*i wil- ..iu ^cc i*£tf ci axicsiBij aavHKittaiii&ia 
j;*.£ lo^u except *or cansaaC taii-cip.*:.. o*: lo^- occurrifijj i.*- a sasw 

Chi generalized version of £ha Ltfcviiio ChaKrfen oayc chat constant 
multiples of special functions HA bo included in the integral as veil aa 
coniUrt: multiples of logs. la aay particular situation aa& yauld be look- 
ins far integrals involving a g-Viii ace o£ cpoci^l £uaccioas only; say, 
eirtor inactions and exponential integrals, is such » case the gaawaiixwi 
version o£ the Liouviila theorem particularizes co give rhe fom of rhe 
in cegr.il as 

q r k 

w»$r$ w, v. , u, . £ad £, ere r&tic^tl egressions in the exponential* log* 

j * 

aritbaia and special fraction cerw already in Che iACdgrznd, e^id che 
e, , d^, cad e. are constants. 

It Id important co note how constant multiples or special funcciens 
arise in the integral. ,Tke sicuizioa ih che log case is quite simple — 
r*-.c intigrind susi he ci che fara 
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1st us recall fran [2] that Che integrand can be expressed ir. cha 



parkin; fraction cccc^posizion 






t*hwre •', in ^«- S^ etc polymzii^li fc - ^ : „_ visSilly £z _. ^ l..-« .\~._ 
with coa-ficieacs involving z&c saititoil, U^fliaa&ia, ^r.w c™j «^uc:«. 
function terns ta the intc-rc^yl. It * c lfc(1; .- -; : — r ^ * ^ c _ ze=^ ift-iclvir.; 

tho uLn v&ri&blc (If any) must &ri*6 firsa wjs*. rczio:;;.! t2r.x. 






Tho other new log c^naa (i.e., not ajnmWlng cj* =^tn variable) nri3o 
through the recursive use oz the Ift^ogtCCiOa process in the poly&Gxial 
pESCi Such terns »;uiC also be obtained frax rctiont.1 :era& at same poise 
i& the recursion. 

2he fort: of the integrand leccinj to She error function zczui is 



I Ice - u 1 » c-^* en(u) 1. 



2as£t on this torn it is clear that the male variable at the point in the 
CatflSrwlOP process when the error ruction arose was an exponential. Further- 
nOM| the tern arose when integrating tha polynomial pert of the asore^cion 
or rational terp* of the forai 
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;?he litzar zeras «e created l£*e polyjwatfcl urn in tha algorithm, whan 
i" is the Min var labia), 

*;^ ah^ll ncu ir<dicct* hoy ^ha aZ^oriths; c^n be axt»£a£ to handle the 
erzer ;3&atiaa« ^st us iw»ue zh^ t::^(i)) is the emit, variable ia the 



kiiUJ^k>Miv j *.»£• — t*L. ... »^_^- »iTi*- a" >fc fc j ^». 



-th. 



i*l s i 



where P, R , find S. are ?oiyaor.i*^ i^ erSfcft;)}* 

For rational terns a reduction ia ti* d&fraa or rfla «^«cr£.ior ccr« be 
;asda. Xhu? for i> 1, 



J?r^-^* |4r*. 
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Uh*r« A *nd S are polynomials in arf;u(«)) abtcinad by c slnpia ccapucation 
based or. the gcd of S find 5*. Par i=l, rtie denominator is completely 
factored zzc we proceed as in J2] ro rind iogfiMcSmie tern*. 
For axosplaj consider tie integral 



C _-kV arffal - i \ 



lie reduction procedure should yia~:d 
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left letter cara need not be factored since the der.ttr.ir.i:tor ia linear 
ir> crf<;:). Ihli tars yields 



-~ io- erf In). 



£u zau L.\Z2„ 



t„M+ .«-* ~ J. 






>:ow leu us consider the polyaesiiil ;*^r^- Sl?p&3ii fcUa oz?£aa&ic& zg at 






vhtre tee A* do aot involve erz(u). Srta the itttegralj i£ it tecUea, mst be 

oi the rora 
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a D+: * r£ (u) * * ft erfn f u > *"** 3 o ■ 



ulcere t£i$ 3- <3o noc involve erf^u), m<1 oaiy 3„ ciey concaixL ne# ioc^ri^h^ic 
or criro;;* function cans. 

We ccn obtain values £o- ;he L. <or decide she nor* in cegrab titty of rhe 

j. 

OJt^TMiioo in logarithmic, exponential, or error-function cema) ay s process 
eiailfil co cho one described in [2]- Au axfin^le which uses chis method Is 
glvcd aelow. Consider 

r -2 

! err (x) dx 

i^criM".; 3. err(x) + 3^ er^C*} -r 3. «r£(x) -r 3- 

« 

£ftMS»adi Aj arz <x) -> A ; erf (^> 4- A- j A- - 1, ^ » 0, A- = & t 



>--.*o,^--L-.*---*-j £^-i g*. * > £ »»i --*♦* *-* *..- ^■•*-i_ l _r^L a ;<J uj ,;*•**". 
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3y matching coefficients, scarciss wico IL, va ob^Lia b ■ « :. Tharafon, 
3* is a cOMCaat, say b^. 



7S b s'"" *V» 



By iatt£»tlAg (trivially on tha leic-siic, by a recursive call to the 
integration procc^ura o& Che righz), wc obcain 



3b* eri (x) +3* c x + co^s££nt. 



So b = 0» 3« = x t b 9 



.2 
(X -r \) A t3/-0 



— U. -r * 2 J C ^ ^ 
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£ -x . T , Ax -x" 

I °2 a + V " * 
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Sag £?_£&) t 3, ■ -77 i t ccai 
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Iczearnrina 



b. er£(x) + L_ = - [2 er£(>:/2) + eoastuc 
1 ° * Ti 



(Saa below for a discussioa of such in iszc^racioa step). 



So b x - 0, B fl o -^3 er£(x /2) + b ( 



Substituting thoso rasulcs into die s^^-*-i £or* of the integral, va 
ObtfclA 

x aui 2 <x) + -^r e" X arifr) - y| er£(x /5) + b Q 

Ic ;hus appears chat whoa the aain vexifible is erf (ufc)) the algorithm 
ij ruito similar to tha log case as we hsv* rlfciaod* The algoriuhn zusc also 
be nncHftafl in the expsnen^iii ca&e, tta have _o pamit constant ciulcipiea 
of uror-<£bCCClOa earns co be included ia the Mluiton. This follows iron 
cur previous diaeuaaion. 2h*t i*, i£ th& E*ir. variable is * * then 



r AMi toW ^ - s^w * a arffec^i, 






Zy ii-iJcrzczia^z;-. vo o,;:«: 



Ao** - (L 1 *10u") ^-ra 
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-^J»l*jy^ -Tf? 1 ) 



Hcr.ca 






We era given A, u, czd k> Our object is » fine L and c, vhsn chaso 

2 
exisc, whera -p - ku. 3 is rcscriczed Co betas rational la die lOsAftthttiCj 

exponential, and error-function rarius which A contains. >icraovar, since we 
do sot tilow algebraic agressions e^ree?: for algebraic constants, we shall 
restrict p accordingly. Sou chat -is csn bo iKpraflafid as a ri*l6nai depres- 
sion in zhe variable of integration zzd iny additional defiae*2 variables 
reprosentias logs, exponentials and error- functions. We want to itnow if -Uu 
is a percict square, This can be decided by factoring. If -ku is a parfact 
rcucre, iar p be a square root o£ i:, Otherwise ignore the error cu&Ctloa 
zarzs and proceed a$ usual. 



Rearranging tha equation abGve, wo obtain 



B 1 + kBu 1 »A • r? 



>J* can now use Che ^achinary of the Sifdh algorithm to solvo for 3 and c, 
i: dwy mist. A solution ror h and c pUflt be cnwnia for if ic vara not it 
rOtrla iCTly zhat an orror-runczion zora is algflbratcelly dependent on an 
ciipononrlal tad zar^s of lover e«?iffl£ty >;a shall assume than exenz far 



:- — *-"■;•' zot^^t.z c»;jc-j:j.^ *«;, -u*, -*-*... *-»> ..*: _ ^^, 



^ in oVG&alBj Itic ^ ccwU-; 



• e -- x 



_ r *-- 



I — - t — £x = La -re ar£(p;. 
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iiere -? * -2j;", so let p = /2 

■ 

DixicrGa elating 



-2x 2 „ ,- t A _ . -2:; 2 _ 2/2 -2x 2 



:- 






Z /it 



L mu*t be tiiiocxl in X* Its dtnoaltmtor is k. The ruir&erator Is 



Substituting 



9o 3 



4-*. «.-!£* 



^0 
£ « -i t 4 - ■ y »^ c. so c = - /2 /^ * * -4i»t * 0, 30 b - 



*J-;j Integral therefore is 



-fi 



-2x 



-/2 A arifr/2) 



ff* ttall nw discuss =:-.* ^rabies* KewafeenA l^ seasrtliiir^ our 



lajoSCiiS *-c£tf- -£* *± ;orw t -u:-u -» — _\_ — + '_ (tta i rf -:.- - _.£ ftu^uar l_ 
iiufiil fi^iiiii^l £uQG£iOai] in Lias;: 2*,fiuW ai » 4 }>L*C£.::£0tt is ifylitivj;.;/ -^__-_ 
cti£ not vary large. 

To* u»h special iur.ee lor, Gna m»jc pa^xars an analysis o£ tac waya in 
which constant xuitiplos of an instance a:T it can encor into the integral. 
This InvolVU some pattern recognition* and in saauis difficult to cbia*f* t 
simp la algorithm to perform the cask in s*^""' Kowaver, £or etch special 
function dsfinsble by an tatafirsl that we hivo uncounrorod, the analysis ia 
quite simple [see 5]* 

He :aust also know in what ways the sp&ciil function interacts uith 
praviou* functions The algorithm assumes the algebraic independence o£ all 
tha tertis. Hence, we mist not let algebraically dependant expressions 
-; r -u: in tn& incegr&nc. Ait c;cr,?:t of such a dependence is given below* 

Jlfc) - 

r 

Ic^zz ra^y *:lso be algebraic relationships between different instances of 
saq scrae function (e.g., log (eh) - log a + log b 4* 2!cnl), In these cat as 
v*i- presently have to rely on previous analyses of the simplification rules 
at£ocintad with the special functions under consideration. Divining the 
.-ii^iiticatioa rules applicable to tny given special function raiacive 
K c givan set of functions taenj; very difficult at present. 

Zz shcls also be clonr that tha fact that an expression is not integrable 
-a isms cf existing functions does nor nscsi-ariLy c*ta tta integral a 
suiCcbjj candidate for a ipsciil function* ?©£ assc?!** noca that 



; d:: 



I.ecause of this relationship zha ia^-rcl m u.;c *oir is c b^d ccndi^tfc 
for a special function. A little* conjidar^ion will revezl sev^r^l restric- 
tions on what oae Pight ceil "primitive" i:\zt£r£i*. S*ficri.c£i;-,& ;pcci.al. 
functions to those defined by &uc;. srt^itivi inscsreis will 3ikc the expan- 
sion e good deal easier. 

L 

In order to uebug the pros~^» wa *»3oed our varcion ax the ai&ch algorithm 
extended to haadie the error function Co chock a fine table of incc^.^U of 
tha irror inaction compiled by \j* 3. Ksurer ia 1953 [6j. (A acre recent 
table, emphasizing definite iniosr^ii; is i?]). Sift program unfortunately 
cannot check the correctness of the entire set of icte^rels because scce 
e^iries involve dlgsbrsic expressions such £S /x end others involve definite 
li^tc^-als. The progrcs: has so fir uncovered a ai&priac in 

f ,2 3 /x 2 i\ -2 , x -x 2 „ 1 -2x 2 

J X «r* dx - (g- - z ) en- x ^^ e erz x -r j^e , 

b&OTA the underlined «x?^c&sion Wta 2 /*:. 

A potentially controversial disagreement between the program znd the 



[ * erf (ox + b) e ?x dx - ^j orf (ax * b) a ?s (2 - 2px + ? 2 x 2 ) 

? 



2 
S- 2 



£."?* /n 



a * arf(aK * b - sr) 



kJt? 



SMe underline*- *oru 4 ; ^.iL «JiL &■.».* -^-- 4 * , ..uVuVjr, ;.r.u:L: ^Lciis*.-*^ 
ir-ciccuQi that ha did cot ahua',t sac 2h-Sf£ u^.v* _-. ^.l *^j^^*i-» ;,v «v;-vc 
decided so: to abaclc Cite integral ait£;ar- a^a £.,*»:, a?d E4irfr.j0M ytlccmt 
to Gn:cr the arena uich a certification of cither of the propo-ed eolation* 
or with another solution of their enooair.^. 

Naurer rec^fctly infonaad u& chac he £panz at LaaoG forty hours ch^c-vlr.^ 
by hand the solution co Che above jro^rfn hafara at save up on it, Efta ?33-6 
at Iirr 3penr ei^ht dilutes en th* problem. He o&lievo that the caachina'i t3oa 
can bo reduced to a tow seconds by caltin;; changes to the sit&plifier usei In 
eh* program- Currently tie siaplifier rare? resents expressions quite frequently 
due to a fear that now logarithmic or error rVncciGn terms were obtained 
in the letter levels of the recursion* Tae algebraic nan ipu le t ions in the 
pro^ra^ ere costly performed using a rational function representation, but 
this representation does not allow one to introduce nsv variables. When a 
new variable appears, the old expressions must be rerepresented to account 
for it& existence. Fortunately the si=^li£iar T G fears are usually ground- 
less, ana appropriate laccifications to it should yield considerable savings 
in execution time. 

Risen has pointed out to us that Ostrouski [3] gives an algorithm which 
-i s La ^hs case where no algebraic expressions appear, en sxteftsioa of the 
T.i&cn algorithm for the log case. Cstro%"£ki caals with integrals or rational 
expressions involving a function defined in terss of an integral. EmavftT, 
he considers only a single function (a.^., log(x} or erf(x •?• 1} ) rather 
than a family of functions. Hence he deals with, the case when a I05 or 
special function is the main variable of integration, but does not handle 
cades vhora constat aulti^lafi of special functi^ Ha appear in the inc££ral. 
jU '.:^ have seen, wis problems IH&C re»ls tfssoivad — such cases ara assan- 
tially sicplii'lcazlcr. jro'atKS. 
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